Problems with Pencils: Lensing Covariance of Supernova Distance Measurements 
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While luminosity distances from Type la supernovae (SNe) provide a powerful probe of cosmo- 
logical parameters, the accuracy with which these distances can be measured is limited by cosmic 
magnification due to gravitational lensing by the intervening large-scale structure. Spatial cluster- 
ing of foreground mass fluctuations leads to correlated errors in distance estimates from SNe. By 
including the full covariance matrix of supernova distance measurements, we show that a future 
survey covering more than a few square degrees on the sky, and assuming a total of ~2000 SNe, will 
be largely unaffected by covariance noise. "Pencil beam" surveys with small fields of view, however, 
will be prone to the lensing covariance, leading to potentially significant degradations in cosmological 
parameter estimates. For a survey with 30 arcmin mean separation between SNe, lensing covariance 
leads to a ~45% increase in the expected errors in dark energy parameters compared to fully ne- 
glecting lensing, and a ~20% increase compared to including just the lensing variance. Given that 
the lensing covariance is cosmology dependent and cannot be mapped out sufficiently accurately 
with direct weak lensing observations, surveys with small mean SN separation must incorporate the 
effects of lensing covariance, including its dependence on the cosmological parameters. 

PACS numbers: PACS number(s): 98.80.-k, 98.62.Py, 98.80.Es, 97.60.Bw 
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Introduction. — Type la supernovae (SNe) have proven 
to be powerful probes of the expansion history of the 
universe 0, contributing to the discovery that this ex- 
pansion is accelerating. A mysterious dark energy com- 
ponent that comprises ~70% of the energy density of the 
universe is presumed to be responsible for this accelera- 
tion. While the presence of dark energy is by now well 
established, its properties and provenance remain a com- 
plete mystery. As the precise nature of the dark energy 
has profound implications for both cosmology and par- 
ticle physics, the elucidation of its properties is one of 
the foremost observational and theoretical challenges. It 
is hoped that more accurate cosmological measurements 
will constrain parameters describing dark energy, and 
eventually shed light on the underlying physical mech- 
anism Several ongoing programs, including the Su- 
pernova Legacy Survey 1 , Carnegie Supernova Project 2 , 
Essence 3 , Sloan Supernova Survey, and Supernova Fac- 
tory 4 , are underway to observe large samples of low, in- 
termediate, and high-redshift SNe and thereby obtain 
~ 10% constraints on the equation of state parameter of 
dark energy. Future attempts to measure crucial prop- 
erties of the dark energy, such as its time evolution, in- 
clude a dedicated space-based instrument as part of the 
NASA/DOE Joint Dark Energy Mission (JDEM). 

It is well known that gravitational lensing provides 
a limit to the accuracy with which the true luminos- 



1 http:/ /www. cfht.hawaii.edu/SNLS/ 

2 http:/ /cspl.lco.cl/~cspuserl/CSP.html 

3 http://www.ctio.noao.edu/~wsne/ 

4 http://snfactory.lbl.gov 



ity distance can be determined for an individual SN 
The effect comes from the slight modification of the ob- 
served SN flux due to lensing by the intervening large- 
scale structure. In fact, the total error budget for SNe 
at redshifts higher than z ~ 1 will have statistical errors 
due to lensing comparable to the intrinsic luminosity dis- 
tance dispersion These lensing effects may have al- 
ready been detected in the current supernova sample [|J , 
although the evidence is still inconclusive |2j. Assuming 
that lensing contributes to the variance of the observed 
SN luminosity distribution (i.e. affects each SN obser- 
vation individually) and using the expected distribution 
function for the cosmic magnification y| , it has been sug- 
gested that the intrinsic power of SNe la observation can 
be restored in the presence of lensing provided the SN 
sample is increased by a factor of 2-3 [5j- 

In addition to the increased variance of SN distance 
measurements due to lensing, spatial fluctuations in the 
foreground mass structures will lead to correlation of dis- 
tance estimates of SNe. Even SNe that are widely sep- 
arated in the radial direction will be lensed by common 
(sufficiently large-scale) modes of the foreground mass 
distribution. 

In principle, one can use fluctuations of the mean in- 
trinsic luminosity to measure magnification statistics [jj . 
While such measurements are useful in the context of 
weak lensing studies, lensing correlations provide a signif- 
icant challenge for precision measurement of dark energy 
properties. The additional covariance due to lensing can 
lead to significant degradation of cosmological parameter 
estimates for future small-field SN searches. It is to be 
emphasized that our results apply to any standard-candle 
approach (e.g., gravitational- wave standard sirens 0]). 
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FIG. 1: Covariance of weak lensing convergence, C K (zi, Zi,9), as a function of two source redshifts, zi and zi, and their 
projected angular separation 9. Left panel: covariance as a function of 9 for several values of z a = z\ — z<x- Right panel: 
covariance as a function of source redshift with the other source fixed at z — 1.7 and for several illustrative values of 9. For 
comparison we also show the lensing variance as a function of redshift. The two horizontal lines represent an intrinsic SN 
measurement error of 0.10 and 0.15 mag (or 0.046 and 0.069 in ScIl/iIl) respectively. Note that lensing variance becomes 
comparable to intrinsic dispersion at z > 1.2 for <7i nt = 0.1 mag, and z > 1.7 for <7i n t = 0.15 mag. The lensing variance at 
low redshift may become smaller than covariance of closely separated SNe when one SN is at high redshift. Of course, the 
correlation coefficient is always less than unity, but can be more than 0.5 if the SNe are separated by ten arcminutes or less. 



Calculational Method. — In order to quantify these 
statements, we first summarize the lensing magnification 
of background SNe due to the foreground mass distri- 
bution and estimate the full covariance matrix associ- 
ated with lensing. Lensing modifies the true SN flux 
by a magnification factor fi, so that the observed flux 
is given by / obs (n, z) = /i(n, z)/ truc (z), where n repre- 
sents the direction of the SN on the sky. In the weak 
lensing limit, this magnification can be related to other 
well-known quantities through |ll| 



M 5 



1 + 2k + 3k 2 + | 7 | 2 + ..., (1) 



where k(<^. 1) is the lensing convergence and \j\ = 
vTi + 72 i s the total lensing shear. Since / obs oc d~[ 2 (z), 
where luminosity distance to a source at a redshift 

of z, fluctuations in \x lead to fluctuations in inferred dis- 
tance so that Sd^/dL = —Sfi/2. Ignoring higher order 
terms (which are suppressed by an order of magnitude 
in the lensing variance one can take /i w 1 + 2k 

and relate SN distance fluctuations due to lensing to the 
convergence along the line-of-sight. Thus, the full covari- 
ance matrix of fractional distance estimates for a sample 
of supernovae is 
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where <7; n t is the intrinsic error that affects each SN dis- 
tance measurement. 



Using the angular cross power spectrum of convergence 
between two different redshifts, computed under the Lim- 
ber approximation 
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Here Jo is the 0th order Bessel function of the first kind. 
In Eq. @, Ti and rj are comoving distances correspond- 
ing to SNe at redshifts and Zj respectively, d A is the 
angular diameter distance, and Pd m (k,r) is the three- 
dimensional power spectrum of dark matter evaluated at 
the distance r; we calculate it using the halo model of 
the large-scale structure mass distribution 0] . 

Eq. (J2J defines the full covariance matrix due to lens- 
ing for supernovae at redshifts Zi and Zj with projected 
angular separation of 6ij on the sky. For reference, the 
previously considered excess variance due to lensing cor- 
responds to diagonal elements of Covy with Zi — Zj and 
6ij = 0. In this limit, 3^(16^) — > 1 in Eq. (0J and one 
recovers the variance, cr 2 (z) = J ldlCg K (z)/2w. 
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FIG. 2: Left panel: Expected errors ont« = const (bottom plot) or w a (with a prior on Qm of 0.01; top plot) as a function of the 
side length of the observed field. The two dashed curves show errors in corresponding parameters when lensing is completely 
ignored, and when only the lensing variance is considered. It is apparent that the lensing covariance contributes to the error 
budget appreciably when the size of the field is < ldeg. We show results for two values of <rg that roughly span the currently 
favored values of the amplitude of mass fluctuations and hence the SN lensing covariance. Right panel: The full expected 
constraints projected into the Q.m — w plane (bottom plot; assuming w = const) and wo — w a plane (top plot; with a prior on 
Qm of 0.01) when as = 0.95 and for the cases of no lensing, lensing variance only, and a few selected survey sizes. We have 
assumed a fixed total number of SNe (iV = 1700) throughout, regardless of the parameter set and survey sky coverage. 



In the left panel of Figure 1 we show the covariance 
C K (zi, Zj,6) as a function of = 9^ (which is assumed 
fixed for the moment) for several values of Zj = Zj, while 
in the right panel we show the covariance as a function of 
z\ with the other redshift fixed at Z2 = 1.7. For reference 
we also plot the variance as a function of redshift z and 
compare it to the intrinsic SN magnitude errors of 0.10 
and 0.15 mag, roughly spanning the error expected in 
upcoming surveys. To estimate the resulting effect on 
cosmological parameter estimates, we compute the Fisher 
information matrix 



~ Op a 



dd L {zj) 
dp/3 



(5) 



If the variance of SN distance measurements alone is con- 
sidered, the Fisher matrix reduces to the familiar form, 
with the factors iV(z i )/(crf nt + of cns ) representing the in- 
verse covariance terms; here N(zi) is the number of SNe 
in the redshift bin centered at Zi and of ens is the vari- 
ance due to lensing. With the full covariance matrix 
considered, this simple form no longer holds. Moreover, 
a full iVtot x iVt t Fisher matrix (and not the redshift- 
binned smaller version) is now required in order to obtain 



the cosmological parameter accuracy estimates; however, 
this is not a novel problem since a correct treatment of 
SN calibration uncertainties similarly requires the full 
iVtot x -Ntot ( or even larger) covariance matrix [l^ . Here 
we implicitly neglect information from the cosmological 
parameter dependence of the covariance matrix; there 
would only be significant information in the covariance if 
the off-diagonal terms were comparable to the diagonal 
ones. 

Discussion. — To estimate cosmological parameter 
measurement errors, we assume a survey with 1700 SNe 
distributed uniformly in redshift out to z — 1.7 (roughly 
following Ref. Q)- To speed up the calculation of the 
1700 x 1700 covariance matrix, we compute it in discrete 
redshift bins, stepping by 0.1 in both Zi and Zj. The co- 
variance also depends on the angular separation of SNe, 
and we distribute the SNe randomly in a square field 
whose side (or total area) we are free to change. The his- 
togram of the angular separations is a smooth bell curve 
that peaks at roughly half the field size. 

Figure 2 summarizes the effect of lensing covariance on 
dark energy measurements from the assumed future SN 
survey. We model the evolution of the dark energy equa- 
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tion of state with redshift as w(a) — wo + (l — a)w a where 
a is the scale factor, and consider measurements of four 
parameters: the matter energy density relative to critical, 
^m, w o, w a, and the nuisance parameter M. that com- 
bines the Hubble constant and absolute SN magnitude 
information. Our fiducial model is standard ACDM with 
Qm = 0.3, w — — 1, and w a — 0. Figure 2 (left panel) 
shows the expected errors onm = const (bottom plot) or 
w a (with a prior on fij^f of 0.01; top plot) as a function 
of the size of the observed field. The two dotted curves 
show errors in corresponding parameters when lensing is 
completely ignored, and when solely the lensing variance 
is considered. It is apparent that the lensing covariance 
contributes to the error budget appreciably when the size 
of the field is < 1 deg. Furthermore, the effects of lens- 
ing covariance depend on the fiducial convergence power 
C K {zi, Zj, 6ij), which in turn is sensitive to the amplitude 
of mass fluctuations as (and, to a lesser extent, other cos- 
mological parameters). Since as is somewhat uncertain 
at present, we show results for two values, og — 0.8 and 
ag = 0.95, that roughly span the currently favored values 
of the amplitude of mass fluctuations in the universe. 

Figure 2 (right panel) shows the full expected con- 
straints projected into the CIm — w plane (bottom plot; 
assuming w — const) and wq — w a plane (top plot; with a 
prior on Qj^f of 0.01) with as = 0.95 and for the cases of 
no lensing, lensing variance only, and a few selected sur- 
vey sizes. Again we see that surveys of less than about 
one square degree will suffer from considerable error due 
to lensing covariance. As the mean separation between 
supernovae is increased, off-diagonal terms in the covari- 
ance matrix decrease, and the resulting effect on cosmo- 
logical parameters is reduced. 

Our results can be understood simply in the limit 
of equal off-diagonal covariance terms. In this case, 
the Fisher matrix estimate of error in parameter p a 
is increased by a factor y/l + (N — l)r 2 relative to 
the case with no off-diagonal terms, where r = 
C K (zi, Zj, e)/^C K (zi, Zi, 0)C K (zj, Zj,0) and N is the to- 
tal number of SNe in the sample. With N ~ 2000 or more 
in upcoming searches, parameter errors will increase by 
a factor of \/2 when r ~ 1/VN ~ 0.02. Furthermore, 
Figure 1 reveals that when z\ > 1 correlations are at the 
percent level when SNe are separated by 9 ~ 10 arcmin. 
Note that in order to accurately estimate the errors on 
dark energy parameters one will need to allow for the de- 
pendence of the covariance matrix on imprecisely known 
cosmological parameters that determine the weak lensing 
convergence power spectrum. Accounting for the uncer- 
tainty in the covariance is even more important since 
galaxy shear maps are not useful to correct for lensing- 
modified SN fluxes 0. 

Conclusions. — We have discussed gravitational lensing 
covariance as an additional source of error for cosmo- 
logical surveys utilizing standard candles. Future su- 
pernova surveys that plan ~ 10-20 deg 2 coverage with 



~ 2000 SNe will be largely unaffected by lensing covari- 
ance. Lensing variance remains an issue, but is reduced 
through increased numbers of SNe (about 50 SNe per red- 
shift bin of width 0.1 are necessary to reduce the lensing 
variance so that it is negligible compared to the system- 
atic floor The cosmological parameter accuracies for 
a survey with a rectangular field that is wide in one direc- 
tion and narrow in another may be compromised, since 
the histogram of the angular distribution of SNe now has 
a peak at an angle of order the narrow side of the survey 
(albeit with a very pronounced tail). We find that true 
"pencil beam" surveys with a sky coverage of < 1 deg 2 
in a single field are subject to significant degradation in 
cosmological parameter accuracies due to lensing covari- 
ance. The consideration of lensing covariance thus argues 
against pencils, and in favor of wider-field surveys. 
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